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Abstract 

We have employed our recently developed Green's function formalism to 
study the dielectric behavior of a model membrane, formed by two periodic in- 
terfaces separating two media of different dielectric constants. The Maxwell's 
equations are converted into a surface integral equation; thus it greatly simpli- 
fies the solutions and yields accurate results for membranes of arbitrary shape. 
The integral equation is solved and dielectric dispersion spectrum is obtained 
for a model corrugated membrane. We report a giant dielectric dispersion as 
the amplitude of corrugation becomes large. 
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I. INTRODUCTION 



It is known that lipid bilayers (abound in living cells membranes) exhibit a ripple phase 
of the bilayer- water interface, in a narrow temperature range IB]. The ripple phase is char- 
acterized by permanent wave-like deformations of the interfaces, the origin of which is still 
a widely debated topic. However, it was observed that a giant dielectric dispersion can 
occur in the radiofrequency range, signified by a large dielectric increment in that frequency 
range [Q. In practice, these membrane systems are under intensive research because they 
are responsible for delivery and retention of drugs [[|. 

We aim to obtain some results on the dielectric dispersion spectrum of corrugated mem- 
branes, by using our recently developed Green's function formalism of periodic interfaces 
Q. In that formalism, we obtained an analytic expression for the Greenian that includes 
the effect of periodicity. In this work, we extend the Green's function formalism to compute 
the local field distribution for a lipid bilayer membrane of arbitrary shape, separating two 
media of different dielectric constants. We will calculate the effective dielectric constant of 
the membrane subject to an ac applied field. 

II. GREEN'S FUNCTION FORMALISM 

The Green's function formalism has been published recently ||. Here we re-iterate the 
formalism to establish notation. We will apply the formalism to a single interface and then 
extend to a bilayer membrane. The electrostatic potential satisfies the Laplace's equation: 

V-[e(r)V0(r)] = -47rp(r), (1) 

with standard boundary conditions on the interface, where p(r) is the free charge density, 
e(r) equals e 2 in the host and €\ in the embedding medium. 

Let V\ and V 2 be the volume of the embedding and host medium, separated by an 
interface S. Denoting 9{r) = 1 if r G V\ and otherwise, leads to an integral equation ||: 
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[l-u0(r)]0(r)=0 o (r) + 



dS' [n'-VG(r, r')]0(r'), 



(2) 



where u = 1 — ei/e2, n' is unit normal to S, G(r,r') = l/|r — r'| and 0o is the solution of 
V 2 0o = -4vrp/e 2 . 

Accordingly, our approach aims to solve a surface integral equation for the potential at 



1. step 1: determine <p(r) for all r e 5 by solving Eq.(@), and then 

2. step 2: obtain 0(r) for all r by using Eq.(@) and the results of step 1. 

In step 1, we encounter a singularity when the integration variable r' approaches the point 
of observation r. To circumvent the problem, we take an infinitesimal volume around r and 
perform the surface integral analytically, we find |J 



where "prime" denotes a restricted integration which excludes r' = r. The (surface) integral 
equation (|3|) can be solved for 0(r 6 S). 



Here we apply the integral equation formalism to a periodic interface. Suppose the 
interface profile depends only on x, described by y = f(x), where f(x) is a periodic function 
of x with period L: f(x + L) = f{x). Without loss of generality, we will let L = 1 in 
subsequent studies. Thus medium 1 occupies the space y < f(x) while medium 2 occupies 
the space y > f(x) separated by the interface at y — f(x). The external field is E and 
0o(r) = — E ■ r is the potential. For a periodic system, </>(r) is a periodic function of the 
lattice vector T. In what follows, we adopt similar treatment as the Korringa, Kohn and 
Rostoker (KKR) method || and rewrite the integral equation as: 



the expense of a two-step solution 0: 




(3) 



III. APPLICATION TO A PERIODIC INTERFACE 
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where the integration is performed within a unit cell. The structure Green's function (Gree- 
nian) is given by |J: 

G(t, r') = £ n' ■ V'G(r, r' + T) = £ ^ ' ( " 7 - ~ ^ ' ( 5 ) 
t m |r — r' — mx| 2 

We were able to evaluate the Greenian analytically [|J: 

fYr „■ r' f(r')) - - % ' ] - Z l^M ( 6 ) 

1 ,y > ,I[ }) cos27r(x-x')-cosh27r(y-/(x')) ' [> 

Eq. (P) is a truly remarkable result - the analytic expression is valid for an arbitrary interface 

profile. If the point of observation (x, y) is located at the interface, the Greenian has a finite 

limit as x' -> i: 

G(x, /(*); f(x')\x' - x) = 1+ ^ ))2 . (7) 

We first solve Eq.([|) for the potential 0(x, f(x)) right at the interface: 

(l - |) /(x)) = -£ /(x) + ^ / dx'G(x, /(x); x', fix'mx', f(x')). (8) 

Then we use Eq.(pJ) to find the potential at any arbitrary point (x, y) ^ S, using the potential 
at the interface. 

(1 - u)4>(x, y) = -E y + — / dx'G(x, y; x', f(x'))(/){x' , f(x')), (9) 



47T 

u 

47T 

for (x, y) G V\ and (x, y) ^ V\ respectively. 



x, y) = -E y + / dx'G(x, y; x', f{x'))^{x' : f(x% (10) 



IV. APPLICATION TO CORRUGATED MEMBRANES 

Here we extend the formalism to a bilayer membrane. Consider two interface profiles 
described by y = ft{x), where £ = 0,1 denote the lower and upper interface profiles re- 
spectively. Again ft(x) is a periodic function of x. Thus medium 1 occupies the space 
fo(x) < y < fi(x) while medium 2 occupies the space y < /o(x) and y > fi(x). For the 
upper (lower) profile, n • y > (n • y < 0), thus the Greenian becomes 



C (r ff T VV f ^ M-^'[f^')sm2n(x-x')-smh2n(f t (x)-Mx'))} nn 

The effective dielectric constant e e of the bilayer membrane satisfies the relation: 

e e £ ^ = J v e{v)E y dV = e 2 E V - ue 2 E y dV, 

where V\ is the volume of the embedded medium. As E = — V0, the volume integration can 
be converted into a surface integration by the Green's theorem. Moreover, for the upper 
(lower) profile, n ■ y > (n • y < 0), thus the effective dielectric constant becomes 

ue 2 Vi n 



e e = e 2 



E V Jo 



- dx [<p(x, fi(x)) - (j)(x, f (x))] 



(12) 



To solve the integral equation, we express the potential at an arbitrary point into a mode 
expansion: 

<f>(x,y) = J2 C ik^j(x)£k(y - ft(x)), (13) 

jk 

where ipj(x) and £k(y) are mode functions. The potential on the interfaces suffices: 

<j>{xj t {x)) = Y,c jk U*)Uv) = J2M^), (14) 

jk j 

where Aj = J2k Cjk£,k{fy- Here we make a few remarks on the choice of the mode functions. 
The choice of the mode function is somewhat arbitrary in theory. In practice, these functions 
should be simple and easy to use. Common choice ranges from extended mode functions 
like the Fourier series expansions to localized mode functions like the step and triangular 
functions [j4|. Substituting the mode expansion Eq.(^) into Eq.([|), the coefficients Ai satisfy 
the matrix equation: 



1-^B-^M 

2/ 4tt . 



A 



-E V, 



(15) 



where 



B 



M,. 



dxipi(x)ipj(x) 



dxdx'ipi(x)G tt '(x, f t (x); x', ft^x'))^^'), 



Vi = j dxipi(x)f t (x). 



(16) 
(17) 
(18) 



It should be remarked that the mode functions need not be orthonormal and the matrix B 
is non-diagonal in general. 

V. NUMERICAL RESULTS 

As a model bilayer membrane, we adopt the interface profiles: 

ft(x) = —a (cos 2nx — 0.1 sin47rx) + (t — 0.5), t = 0, 1, 

where a is the amplitude of corrugation, and the sine function is added to upset the reflection 
symmetry about i = |. The width of the bilayer membrane is thus unity. We adopt the 
step functions for the mode expansions: 

9^{x) = 1, Xi — — < x < Xi + — and Q\ (x) = 0, otherwise, 

where h is the width of the step function. In what follows, we adopt 100 step functions both 
for the lower and upper profiles, equally spaced in the unit interval x G (—1/2,1/2]. The 
integrals Eqs.(16)-(18) can be readily performed. 

To study the dielectric behavior, we apply an ac field at a frequency uj. The embedded 
medium inside the membrane has a complex dielectric constant 

6i = £i + — , 

iuj 

where e% and <j\ are the dielectric constant and conductivity of the embedded medium respec- 
tively, with u being the frequency of the applied field. We adopt the following parameters in 
the calculation: e% = 40, <j\ = 160, while €2 = 80. Also let V\/V = 1. The Maxwell- Wagner 
relaxation time of a planar interface is given by: 

r = ex/ax. 

In Fig.|I], we plot (a) the real and (b) imaginary parts of the complex effective dielectric 
constant e e normalized to e 2 as function of frequency for various amplitude of corrugation 



6 



a ranging from 0.1 to 1.0. As is evident from Fig.|l], there is a giant dielectric dispersion as 
the amplitude of corrugation becomes large (a — > 1.0). 

In summary, we have employed the Green's function formalism to study the dielectric 
behavior of a corrugated membrane. The integral equation is solved and the dielectric 
dispersion spectrum is obtained for a periodic corrugated membrane. We should remark 
that the present formalism can readily be generalized to multi-layers systems as well as 
corrugations in two dimensions @. 
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FIGURES 

FIG. 1. Normalized (a) real and (b) imaginary parts of the complex effective dielectric constant 
of a corrugated membrane plotted as function of the frequency of applied field for various amplitudes 
of corrugation. Note that a large dielectric dispersion occurs for large amplitudes of corrugation. 



9 




10 



